Polaronic hole-trapping in doped BaBiOs 
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The present ah initio study shows that in BaBiOs , Bi''"'' sites can trap two holes from the valence 
band to form Bi^"^ cations. The trapping is accompanied by large local lattice distortions, therefore 
the composite particle consisting of the electronic-hole and the local lattice phonon field forms a po- 
laron. Our study clearly shows that even sp elements can trap carriers at lattice sites, if local lattice 
rela^xations are sufficiently large to screen the localised hole. The derived model describes all relevant 
experimental results, and settles the issue of why hole doped BaBiOs remains semiconducting upon 
moderate hole doping. 

PACS numbers: 71.30.+h,71.38.-k,71.45.Lr,78.20.-e 



The nature of insulator-metal(superconductor) tran- 
sitions in oxides has intrigued scientists for several 
decades dsn. Such transitions can be realised by mod- 
ulating the electrical carrier density using electric field 
effects, chemical doping or pressure [J, [HI . Despite the in- 
tense experimental and theoretical efforts, an atomistic 
interpretation of the experimental data is often exceed- 
ingly difficult, since the available first principles methods 
are too approximate. BaBiOa is a prime example for 
such a phenomenon because it is a charge-ordered insu- 
lator, in which bismuth atoms appear in two different 
oxidation states (3+ and 5+). It undergoes an insulator- 
superconductor transition upon hole dopingj^. Here we 
demonstrate that up to 0.25 holes per formula unit, the 
charge-ordered insulating state prevails, but an insulator 
to metal transition is predicted beyond a hole concentra- 
tion of about 0.3 holes per formula unit. We show that 
the transition is controlled by the coupling between the 




FIG. 1: (a) Monoclinically distorted BiOe octahedra empha- 
sising the tilting {t) and breathing (&) structural instabilities 
around inequivalent Bi'^^ (orange) and Bi''^ (blue) atoms, (b) 
Sketch of the simple cubic (sc) small superlattice comprising 8 
BiOe octahedra. (c) Representation of the large face centred 
cubic (fee) supercell (16 BiOe octahedra) adopted to model 
the polaronic states for x — 0.125 and 0.25, with red arrows in- 
dicating the Bravais lattice vectors. Red atoms indicate Bi^^ 
ions converted into Bi^"^ ions upon hole-doping for x = 0.125 
(large balls) and x = 0.25 (large and small balls). 



additional holes trapping on original Bi'^^ sites and the 
oxygen polarisation fields that surround the BiOg octa- 
hedra. The composite particle, composed of a hole plus 
its accompairying polarisation field (distortion) is known 
as a polaron[7[ 

The primitive cell of BaBiOa can be described as 
Baa^Bi^+Bi^+O^, where Bi^+ and 61^+ cations occur 
in equal parts[l, [3|- Indeed, the two different valencies of 
bismuth are manifested in the measured structural prop- 
erties, from which two different bismuth-oxygen bond 
lengths are derived: the shorter bond length is attributed 
to Bi^+-0 bonds, the larger one to Bi^+-0 bonds [i|. At 
low temperatures, the two Bi species are alternatingly 
ordered in a distorted cubic (monoclinic) structure, in 
which Bi^+ is surrounded by Bi'^"'' neighbours (and vice 
versa), as sketched in Fig. [T]^a) and (b). Upon K dop- 
ing, Bai_j;K2;Bi03 undergoes a phase transition to an 
orthorhombic, still perovskite like semiconducting phase, 
at a composition of x 0.12[l3|. Finally for larger K 
concentrations 0.37 < x < 0.53, the compound becomes 
superconducting [^ with no measurable distortions. Ab 
initio modelling for such a case is problematic, since stan- 
dard semi- local exchange-correlation functionals, such as 
the local density approximation or the generalised gradi- 
ent approximation, predict a too small charge dispropor- 
tionation describing BaBiOs as a metal[ll, 12, ^1, in dis- 
agreement with the experimental observations [lj|. Since 
already the parent material is described incorrectly by 
semi-local functionals, density functional theory (DFT) 
calculations for hole doped BaBiOa are unreliable, and 
until today an atomistic understanding of why BaBiOs 
remains semiconducting upon replacing the divalent Ba 
by monovalent K is missing. 

The drawbacks of semi-local density functionals are 



mainly related to sizeable self-interactions [1 5ll 1 6l | . A the- 



ory that is largely self-interaction free is therefore a pre- 
requisite to model BaBiOs, and for this, hybrid function- 
als mixing a fraction of non-local Hartree-Fock exchange 
(typically 25 %) to the otherwise semi-local exchange are 
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FIG. 2: (Colour online) (a-c) Evolution of the density of states 
(DOS) on Bi^+ (orange) and Bi^+ (blue) atoms with K hole 
doping X for the non-polaronic solution (superlattice with 8 
unit cells). The gray shadows indicate the total DOS. (d) 
Evolution of the DOS of hole-doped BaBiOa simulated by 
removing a fraction of an electron [x) leading to a reduction 
of the gap between the Bi^^ and Bi^^ sub-bands, and the 
simultaneous reduction of the Bi-O bond length difference 
h between Bi^+ and Bi^"*" sites. Symbols A represent the 
estimated values of the breathing distortion b computed from 
the Rice and Wang gap equation [22|. 



a promising choice [17j. In the present work, aU calcu- 
lations are based on the Heyd-Scuseria-Ernzerhof (HSE) 
hybrid functional [3l and were carried out using the Vi- 
enna ab-initio simulation package ( VASP) [l^, [IBl ■ This 
functional can be applied to metallic systems, since the 
long range tail of the Coulomb kernel is screened. 

For pure BaBiOa, HSE correctly predicts a semicon- 
ducting state with an indirect band gap of 0.65 eV [see 
Fig. ID^a)], exhibiting Bi'^+ and Bi^+ sites. For these two 
sites, the Bi-0 bond lengths differ by &=:0.18 A, which 
is in excellent agreement with experiment(6=0.17 A)[^. 
The calculated tilting instability t=11.9° correlates also 
very well with the measured value of 11. 2° Q. The den- 
sity of states shown in Fig. HJa) clearly shows that the 
valence band is dominated by orbitals located predomi- 
nantly on Bi'^^ atoms and the surrounding oxygen atoms, 
whereas the conduction band is dominated by orbitals lo- 
cated close to the Bi^+ atoms. It should be noted that 
we keep the simple formal notation of 3+ and 5+ for the 
Bi-ions, although the actual charge transfer is reduced by 
screening, hybridisation and back-donation as a result of 
the selfconsistent treatment. We will first concentrate 
on the hole doped system without polaronic distortion, 
and then discuss how polaronic distortions modify the 
electronic and lattice structure. 

The non-polaronic state was modelled using a super- 
cell comprising 8 pervoskite unit cells [Fig. [ijb)]. For the 
initial structures — before relaxing the structural param- 
eters — the BiOg octahedra were tilted in accordance to 
the monoclinic structure of Fig. [BJa). In this supcrcell, 
we replaced 2 or 4 of the original 8 Ba atoms by K atoms. 
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FIG. 3: (a) Phase diagram and volume (V^) of Bai-xKxBiOs 
upon doping throughout the insulator-metal transition. The 
computed volumes are compared to neutron powder diffrac- 
tion data[l3|. (b) Schematic structure of the Bi-O plane for 
the insulating charge ordered state of pure BaBiOa {x — 0), 
the bipolaronic phases at x = 0.125 and x = 0.25 and the 
metallic regime at x — 0.5, showing the local oxygen (o) 
breathing environment around Bi"^""" (•), Bi^"*" (•), bipolaron 
Bi3+ Bi^+ (•), and Bi'' '^+ (•) sites. Bi-O bond lengths 
between different sites are displayed. 



corresponding to a hole doping concentration of a; = 0.25 
and X ~ 0.5, respectively. 

For X ~ 0.25, the calculations predict that the two sub- 
stitutional K atoms prefer to be arranged along the [001] 
direction resulting in an orthorhombic supercell after re- 
laxation. The predicted lattice parameters agree well 
with the experimentally observed orthorhombic struc- 
ture at a; = 0.25 (theory: a = 6.11 A, 6 = 6.087 A, 
c = 8.606 A; experiment [l^: a = 6.098 A, b ^ 6.086 A, 
c = 8.584 A). For x = 0.5, K layers parallel to the (100) 
planes are preferred. The final relaxed structure is non- 
centrosymmetric and tetragonal, and the structural pa- 
rameters again agree very well with the experimental val- 
ues (theory: a = b = 4.27, c/a = 2.01, experiment 2 if: 
a = 6«4.26, c/a«2.0). 

The density of states shown in Fig. [D^b) shows that 
the Fermi-lcvel moves into the Bi"^"*" dominated valence 
band, and therefore the system becomes metallic. Fur- 
thermore, as a result of the reduced charge dispropor- 
tionation, the Bi^+ states move upwards in energy until 
the valence and conduction bands overlap [Fig. [D^c)]. At 
x = 0.5 charge disproportionation is no longer sustained 
[Ifc). We also simulated this effect by electronic hole 
doping in the primitive cell (Ba2Bi206) finding quantita- 
tively similar results as for K induced hole doping [Fig. 
2(d)], indicating that chemical effects do not play a major 
role. 
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FIG. 4: (a-c) Evolution of the DOS and (d-f) corresponding bandstructure upon K substitution for the polaronic solution. 
The gray shadows indicate the total density of states, and red curves indicate the bipolaronic states as created by doping, (g) 
Comparison between the theoretical and measured imaginary part of the dielectric function for x—0, 0.125 and 0.25. Peak (A) 
corresponds to excitations from Bi'^"'" into Bi^"*" states, peak (B) corresponds to excitations from Bi'^^ into bipolaronic states 
[red curves in panel (b),(c),(e), and (f)]. The experimental curves (o) are taken from Ref . [2^ (a; — 0) and Ref . [2^ {x = 0.21). 
For the calculated curves, the excitation energies have been scaled by a factor of 0.8 to approximately account for excitonic 
effects otherwise not included in the calculations, (h) Charge density corresponding to the bipolaronic band (red line) of panel 
(e) localised around the BiOe octahedron at the converted Bi'^^ — >Bi^'*' atom. 



Although the above description accounts well for the 
structural properties and the change in volume [see 
Fig. 3(a)], the Fermi-level is located in the valence 
band for any sizable hole concentration, resulting in a 
strong intra-band metallic behaviour, in contradiction 
with experiment 23l. |24| . On the other hand, polarons 
allow the stabilisation of an insulating groundstate. 

The isolated polaron was modelled in an fee like super- 
cell with 16 BaBiOa unit cells [Fig. [ijc)] replacing two 
Ba atoms by K atoms {x = 0.125). The BiOe octahedra 
were again initially tilted as in the monoclinic structure. 
Within this setup the trapping of the two holes on one 
former Bi'^+ lattice site emerges naturally during relax- 
ation. Remarkably, hole trapping and local lattice relax- 
ations are strongly linked. Hence we encounter a truly 
polaronic behaviour. Since a single Bi'^''" cation captures 
two holes, we will refer to this Bi lattice site as "bipo- 
laron" from now on0, [1^. The final relaxed model is 
clearly insulating [Fig. [UJb)] with an electronic struc- 
ture that is similar to the one of pure BaBiOa. The two 
holes are confined in one single unoccupied band below 
the bottom of the conduction band [Fig. 111(e)] . The hole 
trapping is accompanied by a fairly large local relaxation 
and since converted sites have only Bi^+ neighbours, they 
are structurally frustrated and do not order within the 
same local breathing environment as the non-frustrated 
original Bi^"*" atoms, as depicted in detail in FiglHJb). 

The spatial charge distribution of the narrow bipola- 
ronic band is shown in Fig. IDJh). Only about 15 % of 
the charge is localised on the Bi site, the rest is partially 
located on the surrounding O atoms and a significant 



amount of charge is also dispersed throughout the cell 
(not shown). This distinguishes sp systems from d or / 
dominated materials. Fig. ID^h) shows that the band is 
made up by an antibonding linear combination between 
the Bi s orbital and the surrounding O p orbitals (nodal 
planes along Bi-O bonds). Hence the lattice distortion 
accompanying the bipolaron is driven by moving the an- 
tibonding Bi-0 state above the Fermi-level. The bonding 
state is occupied and found in the O subband. This re- 
sults in the observed contraction of the O-Bi bond length 
for the converted Bi site (2.13 A) with respect to the non- 
converted Bi'^^ sites (2.29 A). Attempts to stabilise such 
a solution using more approximate techniques such as 
a one-centre LDA+U treatment failed for this material, 
since the Bi s orbitals are too delocalised. 

Interaction between polarons: stabilising the pola- 
ronic solution at a larger hole concentration {x = 0.25) 
turned out to be rather difficult. To minimise bipolaron- 
bipolaron interactions, one might naively expect the 
bipolarons to arrange in a sc 2 x 2 x 2 super-lattice. But 
relaxation from such a starting structure always turned 
into a metallic state [see Fig. IHb)]. In the fee like su- 
percell with 16 BaBiOa unit cells, however, a brief and 
completely unbiased molecular dynamics run, followed 
by relaxation yields an insulating state. A string of 
Bi'^+ atoms along the face diagonal [101] converts into 
Bi^+ atoms, each trapping two holes [Fig. [ijc)]. Similar 
bipolaronic alignments along all three face diagonals are 
found to be energetically almost degenerate, all leading 
towards a similar insulating state as displayed in Figs, 
m^c) and (f). Compared to the metallic state the insulat- 
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ing solution is about 150 meV/polaron lower in energy. 
The sc arrangement of bipolarons [Fig. [Ijb)] was consis- 
tently metallic for symmetry reasons (lack of some Bragg 
reflections for this arrangement) [2^. As a result of the 
small distance between the bipolarons, we find a rather 
strong electronic interaction between them. The bipo- 
laronic subband [Fig. HJc)] is split into a lower broad 
band and a relatively narrow second subband [Fig. IDJf)]. 
The observed band structure is consistent with weakly 
interacting s like impurity states, with the lower (upper) 
subband corresponding to bonding (anti-bonding) linear 
combinations of bipolaronic states. Although the inter- 
action between bipolarons will partially depend on the 
long range order, the strong interaction beyond x = 0.125 
might well be the origin of the experimentally observed 
monoclinic to orthorhombic phase transition at a; w 0.12. 
The strong interaction ultimately also causes an overlap 
between the bipolaronic band and the valence band re- 
sulting in the metallic non-disproportionated state ob- 
served in our calculations at a; = 0.5 [Fig. Efc)]. A 
proper description of the superconducting state emerg- 
ing at these high doping concentrations will require the 
inclusion of electron-phonon coupling. For the insulating 
state, however, we expect only minor modifications due 
to electron-phonon coupling. 

The comparison between the experimental and the cal- 
culated imaginary part of the dielectric function (62) 
shown in Fig. [3)Jg) confirms that our model captures all 
essential features. In pure BaBiOs a single adsorption 
peak (A) is visible corresponding to the charge-ordered 
excitation between the Bi'^+ and Bi^+ sub-bands. At 
X = 0.125, a second peak (B) emerges related to the ex- 
citation from the valence Bi^"*" band into the bipolaronic 
band. It should be noted that disorder and finite temper- 
ature effects (electron-phonon interaction) will broaden 
the theoretical peaks thus improving the agreement with 
experiment. At x = 0.25, the bipolaronic peak (B) in- 
creases significantly in intensity and shifts towards lower 
energy, whereas peak A loses intensity. Considering that 
our calculations neglect disorder and approximate exci- 
tonic effects, our results compare well with experiment 
and qualitatively explain the emergence of peak B in the 
experimental spectra, its low intensity at low doping con- 
centrations, the downshift and increase in intensity at 
larger doping concentrations (bipolaron-bipolaron inter- 
action), as well as the bleaching of the A mode. 

In summary, we have shown that sp like electronic 
holes can trap at otherwise perfect lattice sites, if the 
lattice is sufficiently flexible and capable of screening the 
local hole. One important observation of the present 
study is that sp mediated polarons are not particularly 
well localised (only 15% of the charge localised at the 
Bi polaron site, see Fig. IH^h)). Resultantly, at higher 
concentrations we find bonding and anti-bonding bipo- 
laronic states and strong ordering tendencies. For sp 
systems the opening of a gap is more akin to a "long- 



range" Peierls distortion than to an essentially "local" 
Jahn- Teller effect. Similar features can be expected for 
many multivalent cations. From a computational point, 
our study demonstrates that hybrid functionals describe 
multivalent cations as well as the resultant polaronic lat- 
tice distortions outstandingly well, and represent a major 
and "not just quantitative" step forward in the modelling 
of oxides. We believe that the successful description of 
the insulating oxide Bai-ajK^^BiOs will spur similar re- 
search in a wide class of complex physical phenomena 
such as trapping of electrons and holes in high-k dielec- 
tric materials or organic semiconductors, or confinement 
at the interface between metallic and insulating oxides. 
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To clarify why the sc arrangement of bipolarons is unsta- 
ble, we inspected the bandstructure of the metallic cubic 
cell [Fig. [2b)], and found three valence bands with max- 
ima at r (empty) and strong and similar 1.5 eV down- 
ward dispersion along FX and FL (states at X and L 
are occupied). Although this supercell is orthorhombic, 
the original cubic symmetry of the perovskite unit cell 
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prevails to a large extend in the bandstructure. A semi- 
conducting state can be expected only if states at F cou- 
ple with those at L and/or X in the 8-fold supercell. At 
X = 0.125 an fee like arrangement fuUfiUs this criteria, 
whereas at 2: = 0.25, the second bipolaron must be placed 
such that the Bragg reflections prevail along space diago- 



nals. This can be realised only by ordering all bipolarons 
in stripes along the face diagonals, as indicated by small 
red circles in Fig. [1] whereas bipolaronic conversion of 
the body centred atom, circled green in Fig. [He), de- 
stroys the Bragg reflection along the space diagonal. 



